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Abstract. We report STAR measurements relating to higher flow harmonics 
including the centrality dependence of two- and four-particle cumulants for harmonics 
1 to 6. Two-particle correlation functions vs. Ar] and A<fi are presented for px and 
^ \ number correlations. We find the power spectra (Fourier Transforms of the correlation 

functions) for central collisions drop quickly for higher harmonics. The Ai] dependence 
of «3{2} 2 and the pr and centrality dependence of t>2 and v% are studied. Trends 
' are conistent with expectations from models including hot-spots in the initial energy 

£h ^ density and an expansion phase. We also present V3 and «2{2} 2 — t>2{4} 2 vs. ^/s NN . 
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Studying the conversion of coordinate space anisotropies into momentum space 
anisotropies has provided important insight into the characteristics of the matter created 
in heavy-ion collisions pp. Elliptic flow (t> 2 = (cos 2(0 — ^rp))) has been studied for 
decades to probe the conversion of the elliptic shape of the initial overlap zone into 
momentum space [2]. In 2007, Mishra et. al. proposed the analysis of v\ for all values 
of n as an analogous measurement to the Power Spectrum extracted from the Cosmic 
Microwave Background Radiation [3]. They argued that density inhomogeneities in the 
initial state would lead to non-zero values for higher harmonics including t> 3 . It was 
subsequently pointed out that information on v\ was to a large extent contained within 
already existing two-particle correlations data jl] , and that v n and v n fluctuations would 
provide a natural explanation for the novel features seen in those correlations, such as 
the ridge like [5] and mach-cone like [6] structures. That the ridge could be related to 
flux-tube like structures in the initial state was already argued by Voloshin in 2003 [7]. 
This was bolstered by hydrodynamic calculations carried out within the NEXSPHERIO 
model which showed that fluctuations in the initial conditions lead to a near-side ridge 
correlation and a mach-cone like structure on the away-side [8]. In 2010, Alver and 
Roland used participant eccentricity (e n> p art ) generalized to arbitrary values of n as 
introduced in Ref. [9] and showed that within the AMPT model, the final momentum 
space anisotropy for t> 3 was proportional to the initial £3 >P art [10] ■ This explained the 
previous observation that the AMPT model produced correlations similar to those seen 
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in the data (albeit with smaller amplitudes) [TT] . Although Mishra et. al. proposed 
measuring v n vs. n to observe superhorizon fluctuations, it was subsequently noted that 
higher harmonics of v n could be washed out by viscous effects and that the shape of v n 
vs. n would provide a tool for studying these effects [12]. The effect of viscosity may 
also make it difficult to observe the effect of the acoustic horizon at higher harmonics. 

In this talk, we present measurements from the STAR collaboration [13] relating 
to two-particle correlations and higher harmonics of v n . Previous interpretations of 
two-particle correlations have focused on jet phenomenology [H] and multi-parameter 
fits adapted from p+p data have been used to succesfully describe heavy-ion data 
without recourse to higher harmonics. Here we compare correlations and their Fourier 
Transforms f^{2} to expectations from the conversion of density inhomogeneities in the 
initial conditions into momentum space anisotropies. 

2. Results 





harmonic n harmonic n 



Figure 1. Top panels: v n {2} 2 for 200 GeV Au+Au collisions in 2.5% centrality 
intervals. Bottom left panel: v n {2} 2 versus n for the 0-2.5% most central data. Bottom 
right panel: a n vs. n from 0-5% most central px correlations with 0.7 < A77 < 2.0. 



In Fig. ED (top) we show v n {2} 2 for n = 1,2,3,4, and 5 for 200 GeV Au+Au 
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collisions in 2.5% centrality intervals using the Q-cumulant method [T5] with unit 
weights. Narrower centrality intervals remove the dependence on weights and bin widths. 
Results are shown for like-sign (LS) and charge-independent (CI) pairs evaluated with 
0.15 < pt < 2.0 and \r]\ < 1. For n = 1, the LS results are negative while the CI results 
are positive. Positive v i{2} 2 indicates particles tend to come out in the same hemisphere 
while negative indicates particles tend to be more back-to-back. Although (vx) is 
expected to be symmetric around mid-rapidity and zero at mid-rapidity, fluctuations 
can lead to deviations from zero from event-to-event and momentum conservation can 
lead to negative values of t>i{2} 2 . From very peripheral to 70% central collisions, the 
n = 2 term is dropping as expected for non-flow correlations. Beyond 70% centrality, 
the effects of elliptic flow overcome the non-flow contributions. The shape of t>2{2} 2 
for 0%-70% central is dominated by the elliptic shape of the initial overlap geometry, 
first rising as multiplicity increases, then decreasing as the overlap region becomes more 
symmetric. For n — 3, t>3{2} 2 in peripheral collisions shows the usual decrease with 
increasing centrality expected for non-flow correlations which are diluted as 1/N as 
the system volume increases (where N is the number of superimposed independent 
correlation sources). But for collisions more central than 70-100%, t>3{2} 2 deviates from 
the typical 1/N behavior and exhibits a centrality dependence that echoes that observed 
for n = 2. This is consistent with the almond shape of the initial overlap geometry 
influencing the n = 3 term as well [T6l [T7] . For higher n, the deviations from 1/N 
become more difficult to discern. We've measured the four-particle cumulant w n {4} 4 as 
well (not shown). Our results for t> n {4} 4 are consistent with zero for all but the n = 2 
term. For the 0-2.5% most central bin, t>2{4} 4 is also consistent with zero. These results 
are consistent with expectations for v n oc £ n , P art and a Gaussian approximation for e x 
and By, where e 2 part = e 2 x + e 2 y [H]. 

The bottom panels of Fig.[T]show v n {2} 2 (left) and a n (right) versus n for 0-2.5% and 
0-5% central collisions respectively. The v n {2} 2 results are taken from the Q-cumulant 
analysis which integrates over all pairs of particles within \r)\ < 1 and 0.15 < < 2.0 
GeV/c. The a n results are derived from a Fourier Transform of the 2-D At] vs. A(f> 

correlation function C = — — — Iri^ 1 — - — PtiVti defined to reflect momentum 

currents [HJ [20j EI] . The Fourier Transform is taken for the region 0.7 < Ar] < 2.0 thus 
excluding short range non-flow correlations such as HBT and p+p like jet fragments. 
The two spectra show the same behavior with the n = 1 components suppressed and 
the remainder of the spectra falling off quickly. The drop with n is well described by a 
Gaussian. The Gaussian width for v n {2} 2 vs. n from all pairs independent of charge is 
o~ n = 1.63 ± 0.02 as expected from the observed A0 width of the near-side peak in the 
two-particle correlations [14] . If a n and t>„{2} 2 are sensitive to the initial eccentricity, a\ 
and v i{2} 2 are suppressed because e\ part ~ (due to the constraint that the eccentricity 
is defined relative to the center-of-mass of the participating nucleons [IE]). Higher 
harmonics are expected to be small due to viscous effects. The shape of a n and t> n {2} 2 
are thought to be sensitive probes of viscous effects in heavy ion collisions. 

Many sources of correlations besides those related to flow can contribute to two- 
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Figure 2. Top panels: t>3{2} 2 vs. Arj for all charged hadrons within two centrality 
intervals in 200 GeV Au+Au collisions. Data are fit with a narrow and a wide 
Gaussian. Bottom left panel: The width of the two Gaussians used to fit v^{2] 2 (Arj) 
vs. centrality. Bottom right panel: the amplitude of the wide Gaussian compared to 
a 1/Ap ar t dependence and a Ap art x £3 part dependence. 



particle correlations and v n {2} 2 . More differential analyses can help to reveal the 
mechanism responsible for the correlations; particularly two-particle correlations vs. 
Ar) and A0 and px — Pt or yx — yx corelations [22J . To look at the higher harmonics 
more differentially we also study t>3{2} 2 vs. the pseudo-rapidity separation between 
the two particles Arj. The top panels of Fig. |2] show examples of t>3{2} 2 vs. Arj from 
200 GeV Au+Au collisions for two centrality intervals, 60-70% central (top left) and 
10-20% central (top right). t>3{2} 2 vs. Arj can be well described by a wide and a narrow 
Gaussian peak. The narrow Gaussian is consistent with expectations from Bose-Einstein 
correlations. By mid-central collisions, t>3{2} 2 is dominated by the wide Gaussian peak. 

In the bottom panels of Fig. El we show the width (left) and amplitude (right) of 
the wide Gaussian for positive, negative, or charge independent particle pairs. The 
width of the narrow Gaussian decreases with centrality but the width of the wide 
Gaussian first increases then slightly decreases with centrality. The wide Gaussian 
reaches a width of = 2 near N paxt = 120 where N paxt is estimated using a 
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Glauber Monte Carlo model [23] with a Woods-Saxon distribution for the nuclear density 
distribution and a 42 mb nucleon-nucleon cross-section. In the bottom right panel, we 
show the amplitude of the wide Gaussian. The amplitude initially drops as l/iVp art 
with increasing centrality as expected for an independent superposition of non-flow 
sources. By N paxt = 50, the amplitude deviates sharply from the l/iVp art trend. We 
compare the centrality dependence of the amplitude to the shape of iV part x e| part . We 
observe that the deviation from the l/N paTt trend is well described by the ansatz that 
f 3 {2} 2 (Ar] = 0) oc iVp ar t x £| part where £§ part is calculated from our Glauber Monte 
Carlo as in Ref [10J. The At] dependence, however, has not been described in detail by 
any model we know of. If t>3{2} 2 is related to the initial eccentricity fluctuations, then 
the reduction of t>3{2} 2 at large Arj would presumably require some decoherence of the 
initial state fluctuations at large rapidity separations [24]. Recent work has found such 
a decoherence effect with a hadron and parton cascade model |25j . 
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Figure 3. Top left: the centrality dependence of vs{2} 2 for Arj > 0.6 as a function of 
^/Np^t■ Top right: «2{2} and ^{2} vs. centrality for two different pr intervals (v n {2} 
analyzed using Q-cumulants with Q vectors reconstructed at r\ > 0.5 and 77 < —0.5). 
Bottom left: v\jv\ vs. centrality for different pt ranges. i<2{2} and t>3{2} vs. pr from 
three different centrality intervals. 

In Fig. [3] (top left) we compare t>3{2} 2 for Ar] > 0.6 to N part x e 2 ipart . The At] > 0.6 
cut excludes the narrow Gaussian shown in Fig. [2] but introduces a dependence on 
the Gaussian width. The like-sign, and charge-independent results are compared to 
IVpart x ^3 iPar t an d found to agree, suggesting v 3 {2} 2 may reflect the initial eccentricity. 
At top right we plot v 2 {2} and v 3 {2} vs. centrality for two pr intervals (see [26]). 
These results are measured using the Q-cumulant method with Q vectors reconstructed 
from different sides of the detector (7/ < —0.5 and rj > 0.5). This builds in a Ai] 
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gap of at least At] > 1 thus eliminating non-flow from HBT correlations and p+p like 
jet fragments but nonflow from highly modified jets could persist. t> 2 first rises as the 
system becomes more dense, then drops as the overlap region becomes more symmetric. 
v 3 shows a similar but weaker centrality dependence. In the bottom right panel, we show 
the ratio of vs{2} 2 /v2{2} 2 from the Q-cumulant analysis with rj gaps vs. centrality for 
six pt intervals. The ratio is smallest in peripheral collisions but then increases with 
centrality, rising abruptly in the central bins. The ratio also increases with with 
the most central results in the highest pt intervals reaching a value greater than unity. 
This is also the px range where 0-1% central two-particle correlations vs. Arj and A0 
show clear evidence of higher flow harmonics before background subtraction [27] . Similar 
behavior was noted before in the ratio of f n {2} 2 /e^ part [28], based on STAR correlations 
data [29] . The rise of v 3 above V2 was predicted in several models but the origin of the 
effect has not yet been fully explained. In the bottom right panel we show t>2 and V3 
vs. pt for three different centrality intervals. V3 increases with p? and shows a weaker 
centrality dependence than t> 2 . 

If v 3 and V2 fluctuations are created from initial density fluctuations and pressure 
gradients, one might expect that they are sensitive to the equation-of-state (EOS) of the 
matter created in the collisions. By scanning y/s NN , one may be able to vary the initial 
energy density so that the EOS is soft and the pressure is smaller thus reducing the 
efficiency of transfer between coordinate and momentum space anisotropies. It remains 
to be seen, however, if this expectation is bourne out in full models. In Fig. H] we show 
V3{2} 2 /el part vs. iV part for v /<^ = 7.7,11.5,39,62.4, and 200 GeV. Above iV part of 
about 50, vz{2} 2 / e\ part rises linearly with iVp art . The rise increases with ■ s Js NN . In the 
top right panel, we show the slope of V3,{2} 2 / e 2 part vs. iVp art as a function of y/s NN . 
The slope shows a generally increasing trend with y/s NN . Based on estimates of the 
Bjorken energy density, the softest point in the EOS should reside near y/s NN = 20 
GeV [30]- Clearly more energy points will be of use in this study. We have also 
studied v 2 {2} 2 — v 2 {4:} 2 ~ 5 2 + 2a 2 2 (see [31]) where 5 2 is the contribution from few- 
particle correlations unrelated to a common expansion plane (non-flow), and a V2 is the 
second moment of the t> 2 distribution which should be sensitive to the initial eccentricity 
fluctuations [32]. u 2 {2} 2 — t> 2 {4} 2 also increases with y/s NN with a possible minimum 
near y/s NN = 20 GeV. Future studies at 19.6 GeV and 27 GeV [33] will help determine 
whether these measurements are sensitive to the QCD EOS. 

3. Summary 

We presented two and four-particle cumulants of v n for n = 1,2,3,4,5, and 6. The 
cumulants are consistent with the ansatz v n oc £ n , P art- STAR also presented two-particle 
correlation functions vs. A77 and A0 in very central collisions (0-1%) for particles 
with 2 < px < 5 GeV/c. Those data exhibit correlation structures consistent with a 
large contribution from t> 3 even before background subtraction. Analyses of the v n {2} 
in central collisions in this pt range (p T > 2 GeV/c) indicate that v 3 {2} 2 > v 2 {2} 2 
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Figure 4. Top left: «3{2} 2 from Q-cumulants with rj < 1 scaled by e 
centrality at 5 beam energies. Top right: The slope of v^{2} 2 / e\ part vs. y/s NN . 
Bottom panels v 2 {2} 2 — «2{4} 2 vs. y/s NN for two centralitics. Bottom panels show 
statistical errors only (smaller than points). 



consistent with expectations of models with fluctuations in the initial eccentricity. To 
better understand the higher harmonics, we examined t>3{2} 2 vs. At] and find that it 
can be fit with a narrow and a wide Gaussian. The narrow Gaussian is consistent with 
Bose-Einstein correlations and narrows with centrality. The wide Gaussian broadens 
to a width of = 2. The amplitude of the wide Gaussian first falls with centrality 
approximately like 1/-/V part until iV part of approximately 50, then it scales as N paTt xel part . 
Based on these observations, it appears data favor the ansatz v \ oc £^ part with an 
efficiency for conversion from £ 2 part into v n {2} 2 that drops off with n like a Gaussian 
for low pt particles but that has a peak at n = 3 for intermediate p? particles. 
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